Cavitation Bubble Dynamics inside Liquid Drops in Microgravity 
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[Almost identical to PRL 97, 094502 (2006)] We studied spark-generated cavitation bubbles in- 
side water drops produced in microgravity. High-speed visualizations disclosed unique effects of the 
spherical and nearly isolated liquid volume. In particular, (1) toroidally collapsing bubbles generate 
two liquid jets escaping from the drop, and the "splash jet" discloses a remarkable broadening. (2) 
Shockwaves induce a strong form of secondary cavitation due to the particular shockwave confine- 
ment. This feature offers a novel way to estimate integral shockwave energies in isolated volumes. 
(3) Bubble lifetimes in drops are shorter than in extended volumes in remarkable agreement with 
herein derived corrective terms for the Rayleigh-Plesset equation. 

PACS numbers: 47.55.dp 



Introduction — Hydrodynamic cavitation is a ma- 
jor source of erosion damage in many industrial sys- 
tems, such as fast ship propellers, cryogenic pumps, 
pipelines, and turbines 1]. Yet, controlled cavitation- 
erosion proves a powerful tool for modern technolo- 
gies like ultrasonic cleaning 0, effective salmonella 
destruction j^, and treatment for kidney stones 0|. 
Such erosion is associated with liquid jets and shock- 
waves emitted by collapsing cavitation bubbles, but 
the relative importance of these two processes remains 
a topic of debate 5]. Fundamental understanding re- 
quires studies of single bubbles in different liquid ge- 
ometries, since bubble dynamics strongly depends on 
nearby surfaces by means of boundary conditions im- 
posed on the surrounding pressure field P, 0. Re- 
cent investigations revealed interesting characteristics 
of bubbles collapsing next to flat 7] and curved 8] 
rigid surfaces or flat free surfaces |9|, |lO| . It is thus a 
promising idea to study bubbles inside spherical drops 
and probe their interaction with closed spherical free 
surfaces. However, for centimeter-sized volumes, such 
geometries are inaccessible in the presence of gravity 
and require a microgravity environment, even though 
gravity plays a negligible direct role for most single 
bubble phenomena. 

In this letter, we present the first experimental study 
of single bubble dynamics inside centimetric quasi- 
spherical water drops produced in microgravity. High- 
speed imaging revealed key-implications of nearly iso- 
lated^ finite and spherical liquid volumes on bubble 
evolution and subsequent phenomena. The results will 
be presented in three steps covering (1) jet dynamics, 
(2) shockwave effects and (3) spherical collapse. 

Experimental Setup (for a detailed description and 
figures see 0) — Microgravity was achieved in 
manned parabolic flights conducted by the Euro- 
pean Space Agency ESA (42nd parabolic flight cam- 



paign). The nominal flight manoeuvres offer sequences 
(i.e. parabolas) of 20 s of weightlessness with a residual 
g-jitter of 0.02-0.05 g at typical frequencies of 1-10 Hz. 
Our experiment and all calibrations were carried out 
on 93 parabolas on three consecutive flight days. 

Each parabola was used to run one fully auto- 
mated experimental cycle, consisting in the computer- 
controlled growth of a centimetric water drop and the 
generation of a single bubble inside the drop. The 
liquid was expelled though an injector tube^ specially 
designed to fix and stabilize the ensuing truncated 
drop (Fig. nj. The bubble was generated through a 
spark-discharge between two thin electrodes immersed 
in the drop. A high-speed CCD-camera (Photron Ul- 
tima APX, up to 120000 frames/s) recorded the fast 
bubble dynamics for a sequence of 11 ms. Three pa- 
rameters were independently adjustable: (1) the initial 
drop size [radius Rd,min = 8-13 mm] by changing the 
expelled water volume, (2) the maximal bubble size 
[radius Rb^max = 2-10 mm] by altering the discharge 
energy between 8 and 1000 mJ, (3) the distance d be- 
tween drop center and bubble center by precisely ad- 
justing the electrodes position. We shall usually de- 
scribe a particular setting in the dimensionless param- 
eter space (q^,£), where a = Rb,max/ Rd,min {relative 
bubble radius)^ e = d/Rd^min {eccentricity). 

Jet Dynamics — As a first result, we overview the 
collapse of eccentrically placed cavitation bubbles in 
the parameter range (0.2 < a < 0.6, 0.3 < e < 0.8), 
recorded at 24 000 frames/s (extracts in Fig. CJ. The 
nearby free surface breaks the spherical symmetry 
leading to a toroidal implosion |^ with a fast micro- 
jet (> 50 m/s) directed perpendicularly away from the 
closest surface element (Fig. HJd). This microjet ac- 
celerates the surrounding volume forming a jet that 
escapes from the right surface in the figure at a re- 
duced velocity (6 m/s). In the meantime, a slower 
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splash escapes in the opposite direction (Fig.QJ:). This 
double-jet picture is consistent with estabhshed studies 
of bubbles in the vicinity of free surfaces 0, , and 
provides the first direct visualization of both bubble- 
induced jets escaping from a steady liquid volume. 




FIG. 1: Two liquid jets generated by the toroidal implosion 
of the a cavitation bubble (a=0.3,s=0.45). Movies at [TT|. 

Close investigation of the splash geometry revealed 
a remarkable and reproducible diameter broadening 
compared to similar jets on ground-based flat free 
surfaces (Fig. I3i,c). In the latter case, the narrow- 
ness is due to a very localized high-pressure peak be- 
tween bubble and free surface 9]. We believe that 
this pressure peak is much broader beneath spherical 
free surfaces, since the distance between bubble bound- 
ary and surface varies more smoothly. Such a pressure 
zone broadening suggests an explanation for a broader 
splash and will be subject to forthcoming simulations. 
Additionally, we discovered that the typical crown sur- 
rounding the splash on flat surfaces (Fig. [^h) is absent 
on the spherical surface (Fig.j^t)- For comparison it is 
enlightening to consider cylindrical free surfaces, which 
are straight and circular in two orthogonal directions. 
One might reason that such surfaces exhibit a hybrid 
behavior between (a) and (c). Indeed, recent studies 
of bubbles inside a falling water stream revealed 
a non-circular crown consisting of two spikes in the 
straight direction of the surface, while no crown seems 
present perpendicular to the figure (Fig. ^p). These 
observations definitely disclose non-trivial features in 
splash dynamics and crown formation and promise an 
interesting avenue for future research. 

Shockwave Effects — Many high-speed visualizations 
contain isolated frames showing a short-lived brilliant 
mist (Fig. This mist lasts for roughly 50 fis and 
is visibly composed of sub-millimetric bubbles (here 
called microbubbles) . Strikingly, these mists exactly 
coincide with the instants of predicted shockwave ra- 
diation. Namely, the spark (or primary) shockwave is 
emitted at the spark generation initiating the bubble 
growth, and collapse (or secondary) Shockwaves are 
radiated at the bubble collapse and subsequent col- 




FIG. 2: Splash on several free surfaces. (On (a), electrodes 
enter from the right, (b) was provided by E. Robert ^13..] .) 

lapses of rebound bubbles Q . Although these shock- 
waves were not directly observable with the present 
setup, their exact synchronization with microbubbles 
clearly discloses a causal relation. We claim that the 
mist of microbubbles is a strong form of shockwave- 
induced secondary cavitation 15]. In other words, mi- 
crobubbles are small cavitation bubbles inside the drop 
volume, arising from the excitation of nuclei (i.e. mi- 
croscopic impurities and dissolved gas) at the passage 
of Shockwaves jT^l That the microbubbles were 

indeed cavities was confirmed by showing that their 
size/lifetime ratio matches theoretical cavity life cy- 
cles. 

The surprisingly high abundance of microbubbles 
compared to faint traces seen in ground experiments 
TE] is reasonably explained by a unique shockwave con- 
finement in the case of isolated liquid volumes. The 
shockwave can be refiected many times on the free 
surface and thereby transfer its whole energy to mi- 
crobubbles via multiple excitations. This full energy 
transformation offers a new way to estimate integral 
shockwave energies. The idea consists in integrating 
the energies of all microbubbles, which can be inferred 
from statistical measures of their densities and maxi- 
mal volumes Vmax via E Vmax{Poo -Pv) .18], where 
Poo is the ambient static pressure and pv is the vapor 
pressure. Applying this method to five concrete sam- 
ples showed that 50-70% of the initial electrical energy 
(200 mJ) is contained in the mist of microbubbles as- 
sociated with the spark shock. It follows that 50-70% 
of the spark energy is radiated in a shockwave, and 
we expect to find the remaining 30-50% in the main 
bubble; a prediction which is consistent with the mea- 
sured maximal main bubble volume corresponding to a 
bubble energy of 70 mJ. The actual precision was lim- 
ited by the optical resolution and non- vanishing con- 
tact surfaces (injector tube, electrodes), and could be 
improved in forthcoming studies. 

Spherical Collapse — We shall now investigate the 
collapse of bubbles centered in drops {e = 0), since 
their spherical symmetry allows a direct comparison 
to bubbles in infinite volumes, governed by Rayleigh- 
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FIG. 3: Example of Secondary cavitation: the primary 
Shockwave radiated between (a) and (b) induces a short- 
hved mist of microbubbles. Movies at 

Plesset's theory Bubbles in drops were found 

to collapse faster in excellent agreement with a the- 
oretical extension derived in the following. Note that 
gravity does not play a direct role in Rayleigh-Plesset's 
theory, nor in our extension, but its absence allowed 
the realization of isolated drops with centimetric di- 
ameters. 

In the present approach we neglect surface ten- 
sion, liquid compressibility, viscosity and mass transfer 
across the bubble boundary. We further assume a uni- 
form and constant temperature, and take the pressure 
inside the bubble as the vapor pressure pv- The ensu- 
ing simplified Rayleigh(-Plesset) equation for a spheri- 
cal bubble with radius Rb{t) in an infinite liquid reads 

—Av 3-2 

= -Rb + RbRb (1) 

P ^ 

where p is the liquid density and = p^ — py with 
Poo being the liquid pressure taken far from the bub- 
ble. Integrating (QJ provided the initial conditions 
Rb{0) = Rb,max and Rb{0) = yields the radius Rb{t) 
of the collapsing bubble. The so-called Rayleigh col- 
lapse time TRayu such that R bjTRqyi) = 0, is given by 
TRayi ^ 0. 915 Rb, max yW^p To comparc differ- 

ent bubble collapses, we present their functions Rb{t) 
in normalized scales where TRayi and Rb,max equal one. 
In these scales, Rb{t) can be shown to be indepen- 
dent of Aj9, p and Rb,max- To assert the validity of 
the above model for the present situation, we carried 
out extensive ground-studies with millimetric bubbles 
in extended water volumes. The data revealed excel- 
lent agreement with the theoretical prediction (Fig. ^ 
dashed line and squares), thus justifying the above as- 
sumptions. 

In microgravity, we studied bubbles centered in 
drops with a relative radius a = 0.50 {Rb,max = 
4.07 mm, Rd rnin — 8.15 mm, cabin pressure Poo = 
80 000 Pa, spark energy = 200 mJ). The bubble col- 
lapses were sampled at 50 000 frames/s with an optical 
resolution of 70 /im. Given this high spatiotemporal 
resolution, the results from 5 experimental runs were 
indistinguishable and ensured reproducibility. The col- 
lapse time was measured as 330±10/is and the evolving 
bubble radius Rbif) was reconstituted from the images 
using a model of optical refraction by a water sphere 



Plotted in normalized scales, the microgravity 
data reveals a remarkable shortening of the collapse 
time (Fig.^ circles). 




Normalized time (= t/Tpayi) 



FIG. 4: Bubble collapse, (dashed line) Rayleigh-Plesset 
theory (Q). (squares) Ground-experiment in extended wa- 
ter volumes, (solid line) Modified theory for spherical 
drops, a=0.5 0. (circles) Microgravity-experiment in 
drops. Errors are given by the size of the data points. 

In order to explain this shortening, we shall now 
derive an equation of motion for cavitation bubbles of 
radius Rb{t) centered in spherical drops of radius Rd{t), 
using the above assumptions. Incompressibility links 
Rb{t) and Rd{t) and conditions the radial velocity r, 
as seen by integrating div v = along a radial axis, 

Rd{t) = {Ri^in + RliW (2) 
f = RbRiy (3) 

For the sake of physical interpretation, we derive the 
equation of motion from energy conservation. The ki- 
netic energy results from integrating the energy density 
^pr^ over the whole liquid volume and the potential 
energy is given by integrating the pressure work pdV 
over the volume traversed by the bubble surface and 
the drop surface. Using Q and 

Ekin = 27rpRlRl{l-X), Epot = ^^Ap (4) 

with X{t) = Rb{t)/Rd{t). Invoking the conservation 
law E}^in-\-Epot = and dividing by AirpRbRl yields 
the equation of motion for bubbles centered in drops, 

= ^Rl + R,R, - 2\Rl - XRbRb + ^RlX^ (5) 

where the term in was obtained by substituting Rd 
using relation Q. This equation comprises the sim- 
plified Rayleigh-Plesset equation (QJ and three novel 
corrective terms, which vanish diS Rd ^ oo. Up to the 
factor AiipRbRl terms in A subtract the change 

of kinetic energy associated with the volume outside 
the drop {Rd{t) < r < oo), which needs to be removed 
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from the infinite medium of Rayleigh's model. The re- 
maining corrective term in accounts for the change 
in kinetic energy due to the variation of the drop's ra- 
dius. This term becomes neghgible for relative bubble 
radii a <^ 1 since X{t) < a\/ t. Integrating Q with 
the initial conditions Rh{^) = Rb,max and Rb{0) = 
gives the collapsing bubble radius Rhit)-, which is in ex- 
cellent agreement with the microgravity data (Fig. ^ 
solid line). This proves that the shortened lifetime is 
entirely due to the finite spherical drop volume, and 
validates the derived equation of motion for bubbles 
centered in spherical drops Q. 

Finally, we apply the validated model to predict the 
collapse times of bubbles with arbitrary relative radii 
a. The energy conservation Ekin + Epot = can 
be solved for 7^5, if Eq is identified with the potential 
energy at maximal bubble radius Rb,max^ 



Rh = 



1 



2Ap fRl^rnax 



X 3p 



f h^max -| \ 



(6) 
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the drop's surface velocity i^c^(t) would approach the 
speed of sound for a significant fraction of the collapse 
phase [use (ISj, Q], thus implying a strong pressure 
drop and preventing further acceleration. 

Conclusion — Microgravity conditions proved essen- 
tial in the study of cavitation bubbles inside centri- 
metric water drops. In particular, (1) we obtained the 
first simultaneous visualizations of both cavity-induced 
liquid jets, and discovered a remarkable splash broad- 
ening and missing crown formation on spherical sur- 
faces. (2) Reflected Shockwaves led to a strong form of 
secondary cavitation and offer an original way for the 
measurement of shockwave energies inside isolated vol- 
umes. (3) Our Rayleigh-Plesset-like model for bubbles 
in drops is in excellent agreement with specific experi- 
ments and allowed a general prediction of bubble col- 
lapse times in drops. Forthcoming microgravity studies 
could fruitfully address the limit of large bubble radii 
(a 1), or focus on direct shockwave visualizations. 



Integrating this equation from Rij = Rb^max to R^ = 
yields an analytical expression for the collapse time. 
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where s substitutes Rb/ Rb.max- The collapse time de- 
pends on a through a collapse factor ^(a) plotted in 
Fig. For infinite liquids (a = 0) we consistently re- 
cover the Rayleigh collapse factor ^(0) = 0.915 jl^ . 
^{a) decreases monotonically with increasing relative 
bubble radius a, and tends to for a oo. Even 
though this limit could not be probed in the micro- 
gravity time available, it seems non-physical for it vi- 
olates the model assumptions. In particular, for a > 6 
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FIG. 5: Collapse factor ^(a) as a function of the relative 
bubble radius a. ^{a) is directly proportional to the col- 
lapse time Tcoiiapse if p, and Rb,max are fixed. 
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